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We investigate compact objects formed by dark matter admixed with ordinary matter made of 
neutron star matter and white dwarf material. We consider non-self annihilating dark matter with 
an equation-of-state given by an interacting Fermi gas. We find new stable solutions, dark compact 
planets, with Earth-like masses and radii from few Km to few hundred Km for weakly interacting 
dark matter which are stabilized by the mutual presence of dark matter and compact star matter. 
For the strongly interacting dark matter case, we obtain dark compact planets with Jupiter-like 
masses and radii of few hundred Km. These objects could be detected by observing exoplanets with 
unusually small radii. Moreover, we find that the recently observed 2 Mq pulsars set limits on the 
amount of dark matter inside neutron stars which is, at most, 10“®Mq. 

PACS numbers: 95.35.-|-d; 97.60.Jd; 26.60.-c 


Astrophysical and cosmological observations reveal 
that most of the mass of the universe is in the form 
of dark matter (DM) 0,0- The nature of dark mat¬ 
ter is still elusive. There are direct methods for de¬ 
tecting DM using particle accelerators ii or analyz¬ 
ing DM scattering off nuclear targets in terrestrial de¬ 
tectors 0, such as CDMS-II Q, CREST-II 0, Co- 
GENT H, DAM^LiBRA 0, LUX [iph SIMPLE [uj, 
XENONlO/100 dl], and SuperCDMS [3. 

Apart from these direct searches, constraints on the 
propert ies of DM can be extracted by studying DM stars 
|l4 Il5l| or the effects of DM on compact objects, such as 
white dwarfs and neutron stars. Neutron stars are the 
densest observable objects in the universe and, thus, a 
privileged laboratory for the exploration of matter un¬ 
der extreme conditions. The structure of neutron stars is 
determined by the equation-of-state (EoS), which is well 
constrained up to normal nuclear saturation density [l6| . 
However, the properties of matter at supranormal densi¬ 
ties are unknown while being fundamental for the deter¬ 
mination of the maximum mass of neutron stars. Recent 


accurate measurements of 2M0 neutron stars [17|, ll8| are 
setting some tension among different determinations of 
EoS and, hence, the structure of neutron stars. 

The possible gravitational collapse of a neutron star 
due to accretion of DM can set bounds on the mass of 
DM candidates Also, constraints on DM can 

be obtained from stars that accrete asymmetric DM dur¬ 
ing their lifetime and then collapse into a white dwarf or 
neutron star, inheriting the accumulated DM 2l|. More¬ 
over, the cooling process of compact objects can be af¬ 
fected by the capture of DM, which subsequently anni¬ 
hilates heating the star 22l426l|. At the same time, self- 


annihilating DM accreted onto neutron stars may affect 
significantly their kinematical properties 27| or provide a 
mechanism to see compact objects with long-lived lumps 
of strangelets 28|- Therefore, it is of high interest to 
analyze the effects of DM on compact stellar objects. 

Recently, neutron stars with non-self annihilating DM 
have emerged as an interesting astrophysical scenario, 
where to analyze the gravitational effects of DM onto 
ordinary matter (OM) under extreme conditions 26, [H- 
33|. In [SOj the two-fluid system of mirror DM and normal 


matter coupled through gravity was considered. An ad¬ 
mixture of degenerate DM with masses of ^1 GeV and 
normal matter was studied in [3l| using a general rela¬ 
tivistic two-fluid formalism. Also, in the approach of 
(30| was followed with various masses in the GeV range 
for DM and different interactions. DM admixed white 


dwarfs were discussed in [3J| assuming an ideal degener¬ 
ate Fermi gas for non-self annihilating DM. 

In this paper we investigate compact objects formed 
by DM admixed with neutron star matter but also with 
white dwarf material. We solve the two-fluid system of 
DM and OM coupled gravitationally while performing, 
for the first time, a stability analysis of the DM admixed 
objects for non-self annihilating DM with a particle mass 
of 100 GeV 


35 


. We find new stable solutions of 
the Tolman-Oppenheimer-Volkoff (TOV) equations with 
planet-like objects of Earth-like masses and radii from 
few Km to few hundred Km for weakly interacting DM. 
For the strongly interacting DM case, we obtain planets 
with Jupiter-like masses and radii of few hundred Km. 
These dark compact planets could be detected via exo- 
planet searches [37| or using gravitational microlensing 
[381441 1. The detection of an unsually small radius im- 
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plying an enormously high mass density would reveal a 
compact planet, much in the same way as the first detec¬ 
tion of white dwarfs. Finally, we explore the DM content 
that neutron stars of 2 Mq El [11 might sustain. 

In order to analyze the dark compact objects formed 
by the admixture of OM with DM, one has to solve the 
TOV equations for a system of two different species that 
interact gravitationally. One proceeds by solving the di¬ 
mensionless coupled TOV equations [il for OM and DM, 
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where p' and p' are the dimensionless pressure and en¬ 
ergy density, respectively, defined as p' = P/m'j and 
p' = p/m'j, with m/ being the mass of the fermion, that 
is the neutron mass and dark matter particle mass, re¬ 
spectively. The physical mass and radius for each species 
are R = {Mp/m'j) r and M = (Mp/m'j) m, respectively, 
where Mp is the Planck mass |42l| . 

For OM, we consider neutron star matter given by the 
equation-of-state EoSI from The EoSs obtained in 


44| are constrained by using input from low-energy nu¬ 


clear physics, the high-density limit from perturbative 
QCD, and observational pulsar data. In particular, the 
EoSI is the most compact one with maximum masses 
of 2M0. We, moreover, map EoSI with an inner and 
outer crust EoS using and [46|, respectively. For 
p < 3.3 X 10^ g/cm^ we use the Harrison-Wheeler EoS 
471. Eor DM, the EoS is taken from for a non¬ 


annihilating massive DM particle of 100 GeV, in line with 
several scenarios for DM, such as asymmetric DM, that 
involve massive particles which cannot annihilate with 
themselves [II, 3-M|. Two cases are studied, weakly 
and strongly interacting DM. The strength of the interac¬ 
tion is measured by the strength parameter, y = mf/mj, 
which is defined as the ratio between the fermion mass 
m/ and the interaction mass scale m/ (see Eqs. (34-35) 
in Ref. 42|). For strong interactions, mj ~ 100 MeV (in¬ 
teraction scale related to the exchange of vector mesons) 
while for weak interactions mj ^ 300 GeV (exchange of 
W and Z bosons). Thus, we take y = 10“^ and y = 10^ 
for weakly and strongly interacting DM, respectively. 

For the determination of the dark compact objects, one 
first has to perform an analysis of the stable configura¬ 
tion for both OM and DM. The stability arguments can 
be found, for example, in [I^, where the stability of the 
different radial modes in a star is analyzed. However, 


FIG. 1: Mass and radius of OM and DM objects as a func¬ 
tion of the pressure of OM. The solid lines indicate the sta¬ 
ble regions for OM and DM, separately. The vertical dashed 
lines delimit the common stable regions for both species. The 
calculation is performed for the weakly interacting DM case, 
with a strength parameter y = 10“^, and the ratio of pres¬ 
sures between DM and OM is equal to 10““*. 


we present here a short description of the stability crite¬ 
ria. We consider small radial perturbations of the equi¬ 
librium configuration, which leads to a Sturm-Liouville 
eigenvalue equation |^. By solving this equation, we 
find that the eigenfrequencies of the different modes form 
a discrete hierarchy with n=0,l,2..., being 

real numbers. A negative value of leads to an ex¬ 
ponential growth of the radial perturbation and collapse 
of the star. The determination of the sign of the mode 
results from the analysis of the mass of the star versus 
the mass density or radius [i^ . The extrema in the mass 
versus mass density (pressure) indicates a change of sign 
of the eigenfrequency associated to a certain mode and, 
hence, a change of stability of the star for a given mass 
density (pressure). The derivative of the radius versus 
the mass density at that mass density determines if the 
change of stability takes place for an even mode (nega¬ 
tive derivate) or for an odd mode (positive derivative). 
Thus, starting at low mass densities where all modes are 
positive, one can perform the stability analysis for higher 
mass densities studying the change of sign of the different 
modes while keeping the hierarchy among them. Only 
when all eigenfrequencies are positive, the star will be 
stable. In this way, one can study the stable regions for 
both OM and DM. 

In Fig. [T] we present this analysis for a ratio of the di¬ 
mensionless central pressures for DM and OM equal to 
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FIG. 2: The total mass of the dark compact object as a func¬ 
tion of the observable radius (radius of OM) for different ratios 
of pressures between DM and OM. Two different strengths 
for the interacting DM are used: y = 10“^ (left column) and 
y — (right column). 


p'dm/p'om — weakly interacting DM case 

{y = 10“^) as a function of the OM central pressure. We 
show with solid lines the stable mass-radius regions for 
OM and DM, separately, while the vertical dashed lines 
delimit the common stable regions for both species (de¬ 
noted with ” stable” legend). We distinguish two common 
stable mass-radius regions, one for OM central pressures 
below 10“®MeVfm~^ and another one for OM central 
pressures above 10“^MeVfm“^. A similar analysis has 
been done for other ratios of pressures and the strongly 
interacting DM case, obtaining similar results. 

Once the stability analysis is performed, we study the 
structure of the obtained dark compact objects, calcu¬ 
lating the total mass of the object as a function of the 
visible radius, that is, the radius of OM, in an equiva¬ 
lent manner as done for neutron stars or white dwarfs. 
Neutron stars and white dwarfs exhibit different mass- 
radius relationships due to their composition and, hence, 
EoS. While neutron stars have typical masses of 1-2 Mq 
and radius of 10 Km, white dwarfs are characterised by 
masses of IMq and radii of few thousand Km. 

In Fig. [2] we present the mass-radius relationships of 
the dark compact objects for different values of the ra¬ 
tio of DM versus OM dimensionless central pressures as 
well as the weakly and strongly interacting DM cases. 
For a small ratio of p'^m/Pom — 10”^; observe 
two stable mass-radius configurations, that correspond 
to the neutron-star and white-dwarf branches of OM. As 
the central pressure of DM increases, the neutron-star 
branch becomes unstable and we are left with the white- 
dwarf branch, with densities for OM below the neutron 
drip line but unconventional masses and radii. Indeed, 


for the weakly interacting DM case and p'dm/p'om — 
10^, we obtain a new type of objects, hereafter named 
dark compact planets (DCPs), with Earth-like masses 
of M =10“"^ — 1O“^M0 and radii from few Km to few 
hundred Km. As for the strongly interacting case, these 
DCPs can reach masses similar to the Jupiter mass, 
ranging from 10“^ — IO^^Mq, and have radii of few 
hundred Km. The DM content of these DCPs ranges 
from M =10“^ — 10“®Mq for Earth-like planets, and 
M =10“^ — 10“®Mq for Jupiter-like ones. Interestingly, 
by choosing a mass of few GeVs for the DM particle, we 
obtain DCPs of similar size but much more massive, thus, 
more compact. Note that by rescaling the EoS for the 
crust with the iron mass, we reproduce similar results for 
our DCPs if we increase the pressure of the dark matter 
relatively to the pressure of ordinary matter. Also, by 
taking the iron mass instead of neutron mass, it will give 
a shift in the maximum mass of white dwarfs of the order 
of the binding energy per baryon relative to the nucleon 
mass, i.e. below 1%. 

At this point, one should consider the different possibil¬ 
ities of having compact stellar objects with DM content, 
that is, by means of the accretion mechanisms of DM 
onto neutron stars and white dwarfs as well as by the 
primordial formation of DM clumps surrounded by OM 
due to free streaming. On one hand, one can estimate the 
DM accreted in neutron stars and white dwarfs by using 
the accretion rate of heavy DM particles onto a neutron 


star 
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and a white dwarf 2]| . A maximum limit for 
the total accreted mass of DM of 10 “^^Mq is obtained 
for neutron stars, whereas in white dwarfs the maximum 
total DM mass accreted is 10 “^^Mq. Thus, the accretion 
mechanism cannot explain the existence of a DCP. Free 
streaming estimates of 1 pc ^ , on the other hand, give 
a lower limit for a total DM mass of 10“^Mq. Interest¬ 
ingly, limits on primordial DM objects from Kepler 531 
and by the Experience pour la Recherche d’Objets Som- 
bres (EROS) collaboration rule out the mass range 
of 2 X 10“® to 1O“^M0 and 0.6 x 10“^ to 15M0, re¬ 
spectively, to constitute the entirety of DM in the Milky 
Way. However, it is still not ruled out that a fraction of 
DM form primordially dark compact objects and accrete 
white dwarf material subsequently, thereby producing a 
DCP that becomes observable with visible light. A de¬ 
tailed analysis of the primordial formation of our DCPs 
will be done elsewhere (55| . 

The search of a DCP could be done in a similar manner 
as for extrasolar planets (exoplanets). Most of them have 
been discovered using radial velocity or transit methods. 
The observation of a DCP using these techniques is pro¬ 
vided by measuring an unusually small upper limit to the 
radius from transiting planets. Gravitational microlens- 
ing, as seen by the MAssive Compact Halo Object (MA¬ 
CHO) project [ 3 ^, the Microlensing Observations in As¬ 
trophysics (MOA) [ 3 ^, the Optical Gravitational Tens¬ 
ing Experiment (OGLE) or EROS M|, could be an 
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FIG. 3: The total mass of the dark compact object as a func¬ 
tion of the mass of DM for weakly interacting DM, y = 10“^. 
The full squares correspond to the neutron-star branch, while 
the full circles indicate the white-dwarf branch. The dotted 
lines match the different symbols given by the actual calcula¬ 
tions for different ratios of pressures between DM and OM. 


alternative tool for detecting a DCP. 

It is also interesting to study the effect of DM on the 
2Mq recent observations of neutron stars 17|, [iSj , in or¬ 
der to determine the amount of DM that such massive 
neutron stars can sustain. In Fig. [31 we display the to¬ 
tal mass of the dark compact object as a function of the 
mass of DM for the strength parameter y = 10“^ and for 
the neutron-star and white-dwarf branches. We observe 
that as the DM content reaches masses beyond lO^^M©, 
a neutron star with 2 Mq can not be obtained. There 
is even a more stringent limit in the DM content coming 
from the reduction of the nominal mass of white dwarf 
of IMq. Similar results for the white-dwarf branch are 
found for the strongly interacting DM case, while the 
maximum mass for neutron star is below 2M0 for DM 
masses > 10“®Mq. 
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